function features = extractOneImage(infile, pvalue, p)
% extractOneImage(infile, pvalue, p)
% infile : input file (JPEG)
% pvalue : select range (for example 0.97 for 97%)
% p : Matrix of input
    %% initialize
    outfile = [infile,'.mat'];
    img = jpeg_read(infile);
    coef = img.coef_arrays{1};
    %% select T
    T = selectRange(coef, pvalue);
    coefarr = mythreshold(coef,T);
    features = extractDCTFeatures(coefarr, p, T);
    %saveFeatures(path, T);
end
